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Abstract—Biometric identification provides high level security 
by identifying individual based on the physiological and 
behavioral characteristics. Iris recognition can be considered as 
one of the most reliable and accurate method of biometric 
technology. The iris is an externally visible and protected organ 
whose unique pattern remains stable throughout adult life. The 
feature extraction plays a very important role in iris recognition 
system. In this paper, an iris recognition system based on 
invariant moment features is proposed. Initially, the iris portion 
is segmented from the original image using Circular Hough 
Transform and Canny Edge Detector. Then the segmented image 
is normalized using Daugman’s Rubber Sheet Model. A set of 
seven invariant moment features are extracted from the 
normalized image after decomposition using Haar wavelet. These 
features are invariant to translation, rotation and scaling. Finally 
the matching process is carried out by comparing the input image 
with the template stored in the database using Euclidean 
distance. The performance of iris recognition is tested on Chinese 
Academy of Science Institute of Automation (CASIA) iris 
dataset. 

Keywords—biometrics; euclidean distance; haar wavelet; 
invarient Moment; iris recognition; 


I. INTRODUCTION 


Iris recognition is a method of biometric authentication 
that uses pattern recognition techniques based on images of 
the irises of an individual's eyes. The authentication is based 
on the physiological or behavioural characteristics of human 
being. Among various biometric techniques such as 
(fingerprint, signature, face, palm-print, palm-vein, iris etc.) 
iris is found to be the most reliable and accurate trait for high 
secured applications. It is because iris is unique, stable, and 
non-invasive[ 1]. 


The iris recognition system consists of an automatic 
segmentation system that is based on the edge detector and is 
able to localize the circular iris and pupil region, occluding 
eyelids, eyelashes and reflections [2]. Features are extracted 
with different feature extraction methods to encode the unique 
pattern of the iris into biometric template. The Hamming 
distance was employed for classification of iris templates and 
two templates were found to match if hamming distance is 
greater than a specific threshold. 
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John Daugman, implemented a working automated iris 
recognition system [3]. The Daugman system is presented in 
[3]. Daugman’s method, based on phase analysis, encodes the 
iris texture pattern into a 256 byte iris code by using some two 
dimensional Gabor filters, and taking the Hamming distance 
[3,4] to match the iris code. 


Ma et al. [5] construct a bank of spatial filters whose 
kernels are suitable for use in iris recognition. They have also 
developed a preliminary Gaussian—Hermite moments-based 
method which uses local intensity variations of the iris [5]. 
They recently proposed an improved method based on 
characterizing key local variations [4]. 


Yu et al. [6] use box-counting method to estimate the 
fractal dimensions of the iris, then classify iris images into 
four categories in accordance with the fractal dimensions. 


Lim et. al. [10] also used the wavelet transform to extract 
features from the iris region. Both the Gabor transform and the 
Haar wavelet are considered as the mother wavelet. From 
multi-dimensionally filtering, a feature vector with 87 
dimensions is computed. Since each dimension has a real 
value ranging from -1.0 to +1.0, the feature vector is sign 
quantized so that any positive value is represented by 1 and 
negative value as 0. This results in a compact biometric 
template consisting of only 87 bits. 


Hu. [18] derived six absolute orthogonal invariants and 
one skew orthogonal invariant based upon _ algebraic 
invariants, which are not only independent of position, size 
and orientation but also independent of parallel projection. 
The moment invariants have been proved to be the adequate 
measures for tracing image patterns regarding the images 
translation continuous functions and noise-free. Moment 
invariants have been extensively applied to image pattern 
recognition, image registration and image reconstruction [9]. 


In this proposed work a set of invariant moment features 
are extracted from the normalized image in wavelet domain. 
In section II, the iris image is denoised and normized. A set of 
invariant moment features are extracted from the normalized 
image in Section II. Section IV provides the experimental 
results of the proposed method. Section V concludes with 
some research perspectives. 


© The Author(s) 2015. Open Access. This article is licensed under a Creative Commons Attribution License 4.0 (CC BY 4.0). 


page 5 


II. IRIS IMAGE PREPROCESSING 


The iris recognition system mainly consists of iris 
localization, normalization, feature extraction, and matching. 
Reprocessing of the iris image plays an important role in iris 
recognition system. 


A. Localization 


Iris localization is a process to isolate the iris region from 
the rest of the acquired image. Iris can be approximated by 
two circles, one for iris/sclera boundary and another for 
iris/pupil boundary. The pupil always lies in a dark area in the 
iris region [11]. 


Canny edge detection is used for generate an edge map. 
Then the circular Hough transform for the iris/sclera boundary 
is performed. From the above techniques all important 
parameters are calculated(radius and centre coordinates for 
both pupil and iris circles) [12]. 


After localizing an iris, the inner and outer boundaries are 
detected. In an eye image, the iris may be partially concealed 
by the upper eyelid, the lower eyelid, or the eyelash. Canny 
edge detector is applied to find the edges of the image and 
circular Hough transform is applied to the edge detected 
image to find the centre and radius of the iris using the 
equation (1). 


(x-a) +(y-by =r 


(1) 
With 
x=a+rcosé (2) 
y=b+rsin0 3) 
Where is the centre of iris,(x ,y) is the iris co-ordinate 


pixel and r is the radius of the iris [12]. 


B. Iris Normalization 

As the size of an iris in a captured image always varies, the 
detected iris is normalized. The normalization is used to 
convert the Cartesian coordinate (x,y) to a polar 
coordinates[4,5]. Figure | shows normalization process. 


The remapping of the iris region from (x,y) Cartesian 


coordinates to the normalized non-concentric polar 
representation is modeled as 
I(x(r,0), y(r,0)) > I(r, 8) (4) 
With 
x(r,0) =(1—r)x, (8) + 1x,(8) (5) 
y(r,0) =U —-r)y,(@) +ry,() (6) 


Where x,(@),y,(@) referred the pupillary boundary 


points and x,(@), y,(@) denoted the limbic boundary points. 


The values of r lies in [0,1] and @ locates in the range of [0, 
2n]. 


I(x,y) is the iris region image, (x,y) are the original 
Cartesian coordinates, (r,0) are the corresponding normalized 
polar coordinates, and are the coordinates of the pupil and iris 
boundaries along the 8 direction [6]. 


Fig. 1 Normalization Process 


C. Iris Enhancement 

The normalized iris image has low contrast and may have 
non uniform brightness caused by the position of light sources. 
All these may affect the subsequent processing in feature 
extraction and matching. In order to obtain a more well- 
distributed texture image, intensity variations across the whole 
image is approximated [10]. 


I. INVARIANT MOMENT FEATURE EXTRACTION IN WAVELET 
DOMAIN 


A. Multilevel 2-D wavelet decomposition 


The wavelet decomposition of the matrix X at level N, 
using the wavelet. Outputs are the decomposition vector C and 
the corresponding bookkeeping matrix S. 


Vector C is organized as 


C=LAW)PHO) | VON) | BQ) |. 


H(N - 1)| VIN - 1)| DIN - 1)... (7) 
| HCL) | Vd) | DC) I. 
where A, ( D, are row vectors such that 
A = approximation coefficients 
H = horizontal detail coefficients 
V= (ical detail coefficients 
D = diagonal detail coefficients 
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Each vector is the vector column- wise storage of a matrix. 
Matrix S is such that 

S(,:) = size of approximation coefficients(N). 

SG,:) = size of detail coefficients(N-i+2) 

fori = 2,3 ...N+1 and S(N+2,:) = size(X)... 


For images, algorithm similar to the one-dimensional case 
is possible for two-dimensional wavelets and scaling functions 
obtained from one-dimensional ones by tensor product. 


This kind of two-dimensional DWT leads to a 
decomposition of approximation coefficients at level j in four 
components the approximation at level j+1, and the details in 
three orientations (horizontal, vertical, and diagonal).The 
following chart describes the basic decomposition step for 
images. 


C (3n+1 sections) 


AT |e [Oy Prainelinay] | 
BS Ey 


§ (n+2-by-2) 


Fig. 2 Wavelet Decomposition 
The Haar wavelet's mother wavelet function w(t) can be 
described as (8). 
1 Ost<1/2, 
yw(t),-1 1/2<1r<1, 


0 otherwise 


(8) 
and its scaling function ¢(t) can be described as (9) 
b(t) ( O<r<l, 
t= 
0. otherwise. 
(9) 


The Haar wavelet operates on data by calculating the sums 
and differences of adjacent elements. The Haar wavelet 
operates first on adjacent horizontal elements and then on 
adjacent vertical elements. The enhanced iris image is 
decomposed using Haar wavelet as in figure 2. 


B. Hu’s Invariant Moment 


For a 2D continuous function f(x,y) the moment of order 
(p + q) is defined as 


M 4 = f [xrytre y)dxdy 


—00 —00 


(10) 


For p, q = 0,1,2,... Adapting this to scalar (greyscale). 


A uniqueness theorem in [15] states that if f(x,y) is 
piecewise continuous and has nonzero values only in a finite 
part of the xy plane, moments of all orders exist, and the 
moment sequence (Mpq) is uniquely determined by f(x,y). 
Conversely, (Mpq)uniquely determines f(x,y). In practice, the 
image is summarized with functions of a few lower order 
moments. 


Central moments are defined as 


Mpg = { | (e- 2)? (yy) Fe ded 


—00 —00 


(11) 


10 10 


and y= 
00 00 


where x= are the components of 


the centroid. 


If f(x, y) is a digital image, then (11) equation becomes, 


Hyg = (x=)? (y= y)" FY) 


(12) 
; Hq 
The normalized central moments, denoted 77, = > 
00 
are defined as 
+ 
Where 7 = = rz +1 for p + q=2, 3... 
Nop = The abscissa variance 
1]o2 = The ordinate variance 
1,, = Covariance of ordinate and abscissa 
7),\.= Distribution intensity towards the right side 
compared to tha faft side of abscissa 
1, = Distribution intensity towards the upper side 


compared to the lower side in ordinate 


1x4) = The abscissa skew intensity 


03 = The ordinate skew intensity 
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Every one of the rotation invariant moments extracts a 
characteristic attribute of the image. For example I, represents 
the moment of inertia along the centroid while I7 extracts skew 
invariant properties which are useful in differentiating 
between images which are mirror reflections of each other 
[16]. 


L, =Mx+Mo2> 

1, = (ho —No) +40 1, 

L=GQig=Siiay + Ct =a) 

Ey = Gg t 3th) FO He) 

1 = 6 =9 ay Fd ig thay 
“SG eias Coy Hae 
px thea) Js 

1, = oo — Mor LOPs0 + 2)” — ay + M03)" ] 
+ 477, 39 +2) 721 + Mo3)> 

LT, = (3721 — 03750 + 2) 
[ot HeY —3Ciy they | 
=Bix =e) 1= Gat Os ee) 
(B76 +%y) =e yey I: 


(13) 


A set of seven invariant moments are extracted from the 
decomposed image[16]. 


C. Matching 


Matching the templates is the last process in the iris 
recognition system. The Euclidian Distance is the most 
common technique used to recognize the iris. The features of 
the sample iris are compared with the template feature vector 
of the irises. The Euclidian Distance is measured using the 
following equation (14). 


(14) 


where the p is the dimension of the feature vector, x; is the 
ith component of the sample feature vector, y; is the ith 
component of the template feature vector. 


ITV. EXPERIMENTAL RESULTS 


The Chinese Academy of Sciences Institute of Automation 
(CASIA) eye image database [10] contains 756 greyscale eye 
images with 108 unique eyes or classes and 7 different images 
of each unique eye. Images from each class are taken from 
two sessions with one month interval between sessions. The 
images were captured especially for iris recognition research 
using specialised digital optics developed by the National 
Laboratory of Pattern Recognition, China. 


100 150 20 250 30 3 100 1 20 20 x0 


(a) Input Image 
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(b) Iris Localization 
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(c) Iris Normalization 
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(d) Iris Enhancement 


Fig. 3 Iris Recognition Process 


The resulted similarity value is compared to the threshold 
level, if the similarity value is greater than or equal to 
threshold level the person is matched (authorized), otherwise 
the person is a mismatch (unauthorized). The Euclidean 
distance is used as a similarity measure for matching purpose, 
maximum accuracy of almost 100 % is obtained compared 
with Gabor filter based method. The results are shown in 
tablel. 


Table 1. Performances of the proposed method 


Methodology Recognition 
- rate 
7 (11) 
Existing Method GABOR, HD 99.25 % 
INVARIANT MOMENT 100 % 
Proposed Method BASED ON 
WAVEDEC,ED 


Fig. 4 shows the overall recognition process of iris image. 
Tests are carried out to confirm that iris recognition perform 
accurately as a biometric for recognition of individuals. 
Experiments are also conducted in order to confirm the 
uniqueness of human iris patterns on CASIA dataset. 
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Fig 4. Recognition Result for User Iris Image 


V. CONCLUSION 


Iris recognition is one of the most reliable methods 
available today in biometrics field. In this paper, a novel 
feature extraction method based on invariant moment features 
in wavelet domain is proposed. The normalized iris image is 
decomposed using Haar wavelet. Since the feature extraction 
plays a very important role in iris recognition system, a set of 
seven invariant moments that are invariant to translation, 
rotation and scaling is extracted. Euclidean distance is used to 
compare the input image with template stored in the database. 
The accuracy achieved by the system is very good. CASIA 
dataset is used for testing the proposed method. 
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